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Abstract 


A surface heating function, defined as the ratio of the time de- 
rivative of the mean annual temperature curve to the surface heat balance, 
is computed from the annual temperature range and heat balance data for 
the North American continent. An annual cycle of the surface heat balance 
is then reconstructed from the surface heating function and the annual temper- 
ature curve, and an annual cycle of evaporative plus turbulent heat loss 
is recomputed from the annual cycles of radiation balance and surface heat 
balance for the continent* The implications of these results for long- 
range weather forecasting are discussed. 



Introduction 


Reliable prediction of monthly mean surface temperature anomalies 

(i. e. , departures from the climatological normals) is one of the major 

goals of long-range weather forecasting. To be of practical use, such 

forecasts should probably depict the spatial distribution of expected anom- 
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alies with a resolution no coarser than 10 - 10 km. However, accurate 

long-range prediction on a scale that small appears rather remote at 
this time. It therefore seems reasonable to accept a more modest objec- 
tive and to investigate the possibility of forecasting monthly mean sur- 
face temperatures averaged over much larger areas, e. g., areas of con- 
tinental scale. While such forecasts can be of little practical value, con- 
tinental scale surface temperature prediction may provide a useful test 
problem for the development of physical methods of long-range weather 
prediction. 

The problem of forecasting the surface ground temperature of a 
continent is formally simpler than that of forecasting the surface air tem- 
perature, as no advective terms appear explicitly in the thermodynamic 
equation for the ground. In. view of the crude state of the art, it will there- 
fore be assumed that the ground and surface air are in thermal equilibrium 
on a monthly time scale, and that the two temperatures are identical. 

The physical law governing the variation of surface temperature 
with time, t, may be written in simplest form as 


8u 

at 


= kh, 


where u is the daily surface temperature, h is the daily surface heat 
balance, and k is a corresponding surface heating function. If h is ex- 
pressed in langleys per unit time, k is in units of degrees per langley 
(deg ly *)• The surface heating function, k, may 1 e interpreted as the re- 
ciprocal of the product of an average heat capacity and an "equivalent 
depth of thermal influence" of the continental surface layer. (See Appendix). 


In place of the "daily” temperature one may substitute the tem- 
perature observed on any observational schedule, c. g. , hourly, twice- 
daily, etc. 
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However, it is probably simpler merely to consider equation (1) as 
the definition of k. The surface heat balance, h, may be written as 

h = r - e - p (2) 

where r is the daily surface radiation balance, e the evaporative heat 
loss by the surface, and p the turbulent ("sensible'') heat loss. (See, 
e.g., Sellers, 1965.) 

Let a bar denote the monthly mean value, and a prime the daily 
deviation from the monthly mean. Also, let a capital letter with an 
asterisk denote a climatological monthly normal, while a capital letter 
without asterisk denotes the monthly anomaly. Then, 


u = u + u 1 , 

* 

u =U+U , 
etc. 


h = k + k 1 , 
k = K + K* 


h =R +h*, 
h = H + H*, 


From (1) and (3) 

3}S 

ATT 9U * * * * 

+ r-— = K H +KH +K H+KH+(k'h'). 
ot at 


(3) 


(4) 


It will be assumed that, for the climatological normal values. 


aiT 
TT = 


k" h"' , 


(5) 


which may, in fact, be considered to be a definition of K . Furthermore, 
it is probably safe to assume that (k'h 1 ) may be neglected in (4) and that 
K << K*. Thus, an approximate prognostic equation for the monthly surface 
temperature anomaly is 


atj * * 

— =KH +K H , 


(6) 



which might perhaps provide a basis for the prediction of monthly 
surface tempe r ature anomalies by one-step time extrapolation. 
(Experimental one-step predictions of monthly temperature anomalies 
based on a more elaborate physical model* have been described by 
Adem (1965, 19/0).). 

The evaluation of (6) requires a knowledge of the monthly 
anomalies of the surface heat balance and surface heating function, 

as well as their monthly climatological values. However, an effort 

* * 

(described below) to evaluate K , as well as H , for each month 

of the year from climatological data was not successful, and Indeed 

* 

It proved necessary to assume a constant value of K all year round 

* 

In order to derive a consistent set of monthly H values. In view 
of the need for such crude approximations, it appears futile, at 
this time, to attempt to evaluate K. The approximate prognostic 
equation is thus reduced to 



H. 


* 

where K is treated as a constant all year round. 

The practical applicationsof (7) # including the evaluation 

of H and the extrapolation of U, are beyond the scope of this paper. 

* 

We have, however, attempted to evaluate monthly values of H 
and a mean annual value of K over a continental area, namely 
the greater part of North America. The results of that effort are 
described in this note. 



Calculations and Results 


Maps of mean monthly surface (air) temperature are available 
in a variety of climatological atlases. For the continental mean monthly 
surface temperature calculations we have used the Northern Hemisphere 
atlas of Crutcher and Meserve (1970). The most comprehensive source of 
global surface heat balance data is the atlas of Budyko (1963), in which are 
presented global maps of the mean monthly components of the surface heat 
balance, including the surface radiation balance and the evaporative and turbul 
ent heat losses. 

Mean monthly values of the temperatures and surface heat balance 
components were numerically integrated with respect to area over North 
America (based on interpolations between isopleths and summations over a 
5 degree latitude - longitude grid) between latitudes 30 N and 60N. Budyko 
(1963) has computed the average annual values of the heat balance components 
over all continental and oceanic areas, but does not present the monthly values 

The continental average monthly values of the surface radiation 

£ 5jc Jjc 

balance, R , the evaporative heat loss, E , and the turbulent heat loss, P , 
were computed separately, and the surface heat balance, H , was calculated 
from 


* * 
H = R 


- (E 


+ P ) 


The monthly heat balances were then added to obtain the annual heat balance, 
which should be zero. As shown in Table 1, this sum was not quite zero. 

A monthly correction was therefore computed by dividing the sum by twelve, 
and this constant factor was then subtracted from each monthly value to ob- 
tain a set of adjusted monthly heat balances for which the annual sum is ef- 
fectively zero (except for roundoff). The correction is undoubtedly smaller 
than the errors of estimation, such as those resulting from interpolation be- 
tween isopleths, and docs not exceed 10% for any month. The computed 
values of R , E + P , H , and the adjusted value, H , together with the 

,;c ^ # 

monthly values of U , are shown in Table 1. The annual values of R and 
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** + P in Table 1 are about 10 percent higher than those computed 

by Budyko (19&3) for North America, probably due to our use of 60° N 
as the northern limit of integration. 

The original intention of this study was to compute K for each 

$ 

month from equation (5), with H calculated from (8) and the temperature 

derivative computed from the annual temperature curve; but this was not 

possible. On physical grounds, K should be positive and finite for all 

months. Therefore, H should be zero at the times of maximum and 

minimum monthly mean temperature, with positive maximum during the 

spring warming and negative minimum during the fall cooling, i. e. , 

approximately one season out of phase with the temperature curve. How- 

* * 

ever, as shown in Table 1, the pnase relations between H and 9U / 9t 
appear to be unrealistic, and it is impossible to compute a set of reason- 

51 ' 

able monthly values of K from equation (5) and the tabulated data. For 

* 

example, a direct computation of K would give values of zero twice a 

year, at the times of minimum and maximum temperature, and infinity 

* 

twice a year, when H reverses sign. Furthermore, negative values 

of K would be computed for some months. Part of the difficulty un- 

* 

doubtedly arises because H is computed as the relatively small difference, 

* * v ' 

in most months, between two quantities, R and E * 7 , which are, 

in most months, of the same sign and approximately equal magnitude. 

# 

The computed annual cycle of R appears to be realistic, but that of 
# £ 

E +P is far less reliable. Thus, the difference, H , exhibits an un- 
realistic annual cycle. 

Sj i 

To calculate a constant annual value of K , it was assumed that the 
integral of the computed surface energy balance over the 6 -month period 
from the January temperature minimum to the July temperature maxi- 
mum is a reliable quantity. K was thus assumed to be equal to the ratio 
of the ai nual temperature range, in degrees Celsius, divided by the 
total surface heat balance, in kilolangleys (kly), over the same half-year 
period. (Since the temperature range is in effect calculated between 
mid-months, the half-year heat balance is computed as the sum of the 



monthly values for February through June plus half the sum of the 
values for January and July. ) From Table 1, the mean annual tem- 
perature range from mid-January to mid-July over North America is 
0 

28.61 C, and th<~ half-year surface heat balance is 2.39 kly. Thus, 

* - 1 

K is found to be 12 deg C kly . (The same numerical value is, of 

course, obtained if the calculation is done for the 6-month cooling 

♦ 

period from July to January, as the annual sum of H is effectively zero*) 

Table 1* Monthly computed values of surface radiation balance (R ), 

* # 

evaporative plus turbulent heat loss (E + P ), unadjusted 

♦ # 
surface heat balance, H , and adjusted heat balance, H , in 

-1 a 

kilolangleys per month (kiy mo ), for North America between 

latitudes 30 N and 60 N, derived from Budyko (1963). Also 

shown are monthly mean temperatures, U , in degrees Celsius 
0 

( C), derived from Crutcher and Meserve (1970), for the 
same region. 


Month 

* 

R 

« * 

(E + P ) 

* 

H 

H 

a 

U 

* 


Jan 

- 0. 28 

+ 0. 72 

- 1.00 


0. 96 


9.27 


Feb 

+ 0. 73 

+ 1.49 

- 0. 76 

_ 

0. 72 

• 

7. 57 


Mar 

+ 2. 60 

+ 2. 75 

- 0. 15 

- 

0. 11 

- 

1.45 


Apr 

+ 4. 75 

+ 3. 72 

+ 1. 03 

+ 

1.08 

+ 

5. 47 


May 

+ 6. 74 

+ 5. 73 

+ 1. 01 


l. 06 

+ 

11. 55 


Jun 

+ 8. 07 

+ 6. 86 

+ 1. 21 

+ 

1. 26 

+ 

16. 04 


Jul 

+ 7. 90 

+ 7. 35 

+ 0. 55 

+ 0. 60 

+ 

19. 34 


Aug 

+ 6. 45 

+ 5. 85 

+ 0. 60 

+ 0. 65 

+ 

18. 80 



+ 4. 29 

+ 4. 69 

- 0.40 

• 

0. 36 

+ 

13. 82 


Oct 

+ 2. 19 

+ 3. 07 

- 0. 88 

- 

0. 84 

+ 

8. 21 


Nov 

+ 0. 63 

+ 1.43 

- 0. 80 

- 

0. 76 

+ 

0. 21 


Dec 

- 0. 04 

+ 0. 91 

- J. 95 

• 

0.91 

- 

6. 03 


Sum 

44.03 

44. 57 

- 0. 54 


0. 01 





Correction: 0. 54/12 = 0, 045 











On the assumption that K is constant throughout the year, it 

* 

is now possible to calculate new monthly values of H from equation (5) 

and the monthly mean temperatures. For this calculation, the re- 

♦ 

suits of which are shown in Table 2, 9U /9t was evaluated for each 

month by centered differences over two months. From the new surface 

heat balance, H , and the presumably reliable radiation balances in 
c * * 
Table 1, it is now possible to calculate a new set of values of E + P . 

» * 

These estimates, designated as (E + P ) , are also listed in Table 2. 

♦ # ^ 

The annual cycle of (E + P ) in Table 2 appears to be at least as 

* * c 

reasonable as that of (E + P ) in Table 1. 

The results of the calculations are illustrated in Figure 1, 

* * * * 

showing the annual cycles of U , R , H , and H , and in Figure 2, 

a * * c * * 

in which are plotted the annual cycles of R , H , and (E + P ) . It 

* * c c 

should be noted that the scales for R and H are not the same in Fig- 
ure 1, whereas a uniform scale is used in Figure 2. 

* 

From the new values ox H c in Table 2, the half-year surface 

heat balance is found to be 2. 20 kly rather than 2. 39 kly. If this value 

* -1 

is used to calculate K , the heating function is found to be 13 deg kly 

rather than 12 deg kly~*. 


Discussion 

Equation (7) represents, in its simplest form, the basis for 

thermodynamic prediction of monthly surface temperature anomalies 

* 

over a large continental area. If the estimated value of K (13 deg C 
kly"*) is correct, an anomaly of about ±0. 1 kly mo * in the surface 
heat balance is required to produce a surface tcmpc rature anomaly 
tendency of ± 1 deg C mo“*. This energy flux represents at least 
20 percent of the normal monthly surface heat balance over North 
America. However, a heat balance anomaly of 0. 1 kly mo also * e- 
presents less than 5 percent of the normal monthly surface radiat'on 
balance (or the normal monthly evaporative plus turbulent heat loss) 
over the continent in any of the months from March through October. 



Thus, a forecasting system based on the monitoring of radiative 
fluxes requires very high precision in the radiation measurements. 
An even greater difficulty lies in the evaluation of the evaporative 
plus turbulent heat losses. 

It is apparent from Figure 2 that the monthly normal values 

* * * 

of R and (E + P ) are highly correlated and just slightly out of 
phase. The annual temperature cycle is seen to be the result of only 
a slight imbalance between surface heat sources and sinks. At the 
present time no program exists for monitoring these surface energy 
fluxes on a large scale. Even meteorological satellites are, at best, 
capable of measuring only the radiation balances at satellite altitude. 
Thus, the direct evaluation of monthly anomalies of surface heat bal- 
ance over a continent is clearly a formidable problem, and does not 
appear to represent a viable approach to routine long-range tem- 
perature prediction. However, the thermodynamic method may be 
useful when anomalies in the surface energy fluxes constitute large 
fractions of the normal surface fluxes. Such extreme anomalous 
events are most likely to occur in winter. Unfortunately, very little 
is known about month- to -month and year-to-year variations in the 
surface energy fluxes over continents. Thus, there is as yet little 
or no empirical basis for thermodynamic prediction by the method 
described above. 
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Table 2. Monthly surface heat balance* H * computed from estimated 

* - 1 c 
K (12 deg kly ) and mean annual temperature curve for 

North America. Also shown are monthly estimates of 
* * * * * * 

E + P computed from (E + P ) = R - H >s are 

C s 

kly mo ~ , 


* 

Month ** c 


(E* + P*) c 


Jan 


0.06 

.. 

0. 22 

Feb 

+ 

0. 33 

+ 

0.40 

Mar 

+ 

0.54 

+ 

2. 06 

Apr 

+ 

0. 54 

+ 

4. 21 

May 

+ 

0.44 

+ 

6. 30 

Jun 

+ 

0. 32 

+ 

7. 75 

Jul 

+ 

0. 12 

+ 

7. 78 

Aug 

- 

0. 23 

+ 

O' 

• 

O' 

00 

Sep 

- 

0.*4 

+ 

4. 73 

Oct 

- 

0. 57 

+ 

2. 76 

Nov 

- 

0.59 

+ 

1. 22 

Dec 

- 

0.40 

+ 

0. 36 


Sum 


0.00 


44. 03 
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Figures 


Fig. 1 Annual cycles of mean monthly temperature, H ( C), surface 

* - 1 * 

radiation balance, R (kly mo ), surface heat balance, H 

- 1 ^ 
(kly mo ), derived from Budyko (1963) atlas only, and re- 

* - 1 

computed surface heat balance, H c (kly mo ) derived from 
Budyko data and annual temperature cycle for the North 
America continent between latitudes 30N and 60N. 

Fig. 2 Annual cy os of surface radiation balance, R (kly mo *), 

recomputed surface heat balance, H (kly mo ), and re- 

c * * 

computed evaporative plus turbulent heat loss, (E + P ) c 
for the North American continent between latitudes 30N 
and 60N. 



Appendix 


The heat conduction equation for the surface layer of a continent 
may be written as 


9T 1 ‘ 9£ 
9t ~ C dz 


where T, the earth temperature, and q, the vertical heat flux, are both 
functions of depth, z,as well as time, t, and C is the heat capacity in the 
layer. The heat equation may be integrated from the earth' s surface, 
where q = h, down to a depth, L, where q = 0. Thus, for a homogeneous 
layer, 



o 


]_ 

C 


r L do 

.) 

o 


dz 


h 

C 


For the mean value theorem, 



9T 

9t 


dz 


( 

\ at ! 


L 


where the bar denotes a mean value with respect to depth. Let an 
"equivalent depth of thermal influence", D, be defined such that 



where u denotes the surface temperature. Then, 


9u _ _h_ 
9t ” CD 


kh • 


In an isothermal mixed layer such as may be found in the ocean, 

L = D, while for a continental surface layer, L is greater than D. For 
example, if the amplitude of the temperat ure variation were a linear function 
of depth, vanishing at the depth L, then L would be equal to 2D. However, 



in the more realistic case of an exponentially damped temperature 
variation (which follows from the solution of the heat conduction equation 
for a layer of uniform thermal diffusivity), L is greater than 2D. 

From the annual cycles of temperature and surface heat balance 
averaged over the North American continent, a mean annual heating 
function is computed in this paper to be about 12 degrees Celsius per 

_ 3 

kilolai.gley. The average heat capacity of soils is roughly 1/2 cal cm 
degC”* (Sellers, 1965). Thus, the mean value of D is approximately 
1 1/2 meters and L is greater than 3 meters for the annual cycle. Sellers 
(1965) indicates that the annual temperature cycle over land penetrates 
to depths of 5-20 meters. 
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